Coexistence of diffusive and ballistic transport in a simple spin ladder 
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We show that in a nonintegrable spin ladder system with the XX type of coupling along the legs 
and the XXZ type along the rungs there are invariant subspaces that support balhstic magnetization 
transport. In the complementary subspace the transport is found to be diffusive. This shows that 
(i) quantum chaotic systems can possess ballistic subspaces, and (ii) diffusive and ballistic transport 
modes can coexist in a rather simple nonintegrable model. In the limit of an infinite anisotropy in 
rungs the system studied is equivalent to the one-dimensional Hubbard model. 



Transport properties of simple spin ladders are inter- 
esting for two reasons. On one hand they are realized in 
a number of materials [H-Qj on the other they serve as 
model systems on which theoretical ideas can be tested. 
For instance, the Hubbard model - a paradigmatic model 
of strongly correlated electrons - is equivalent to a lad- 
der system. One of the most actively investigated areas of 
statistical physics in recent years is nonequilibrium prop- 
erties of strongly correlated systems. In particular, there 
is a long quest to understand transport properties of sys- 
tems from first principles, e.g. M. For recent studies of 
transport in spin ladders see |6l4ll|. With rapidly pro- 
gressing experimental cold atoms techniques [12[, trans- 
port properties of Fermi gases as well as of the Hub- 
bard model [31 have actually been measured. Perhaps 
the simplest question, is a given model diffusive or ballis- 
tic, seems to be for most models too difficult to rigorously 
answer, even if the system is integrable, an example be- 
ing for instance the gapped one-dimensional Heisenberg 
model, see, e.g. |5j. A powerful method to prove ballis- 
tic transport is by bounding the time-averaged current 
autocorrelation function using constants of motion, the 



sesses a number of ballistic invariant subspaces. That is, 
there exist subspaces of the total Hilbert space that are 
invariant under unitary evolution generated by H, and 
in which inhomogeneities spread out with a constant ve- 
locity, implying ballistic transport. In addition, we shall 
show that in the rest of the Hilbert space the transport 
is diffusive. The Hamiltonian is, 
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so-called Mazur's inequality 



16|. Because quantum 



The model studied, Eq. ([T]), contains the integrable Hub- 
bard model as a special limit when J — > and JA -> 
const. [l9l |. Let us first briefly mention standard sym- 
metries of the XX ladder. Total magnetization along 
the z-axis, Z = Y^^=i ^| + conserved. There are 

also two lattice symmetries, namely, a parity that 
exchanges sites j and L + 1 — j , and a parity Py that 
exchanges the two chains. In the sector with Z = there 
is an additional spin-fiip symmetry corresponding to the 



integrabihty is usually defined^Jby the existence of an transformation F = n, ^^rf. There exist though addi- 



infinite set of local conserved quantities it often leads to 
ballistic transport. In fact, all proved ballistic systems 
are integrable and possess either local [l6| or quasilo- 
cal [l8| conserved quantities that have nonzero overlap 
with the current. Based on this one is tempted to con- 
clude that ballistic transport is possible only in integrable 
systems and that chaotic ones display diffusion. However, 
as we shall show, this widely held belief is not correct. 
Studying a class of spin ladder systems, which includes 
nonintegrable as well as integrable instances (the Hub- 
bard model) , we explicitly show the existence of ballistic 
transport. This provides a new mechanism of ballistic 
transport, different from the so- far known ballistic trans- 
port in integrable systems. 

XX-ladder. - We shall study the so-called XX spin lad- 
der composed of two coupled spin- 1/2 chains in which 
the nearest-neighbor coupling along two chains (legs) is 
of the XX type, while the interchain coupling (rungs) is 
of the XXZ type. We shall show that the XX ladder, 
regardless of the value of two parameters J and A, pos- 



tional symmetries that are a consequence of an XX-type 
interaction along the legs. Discussion is simple in the 
rung eigenbasis. On one rung the eigenvectors and eigen- 
values are the singlet |S) = (|01) - |10))/%/2 with £'s = 
-(2-1- A) J, and three triplet states, |T) ee (|01) + |10))/%/2 
with St = (2 - A) J, |0) = |00) with Eq = JA and 
|I) = 1 11) with El = JA (here the 1st or 1 in the 
ket/bra denotes the state on the upper leg, the 2nd on 
the lower leg). Written in the rung eigenbasis, labeled 
by 4 letters {S,T, 0,1}, the leg nearest-neighbor Hamil- 
tonian has two kinds of terms. One that represent 
hopping, and the other that cause "scattering" . The hop- 
ping terms are 

|0S) o ISO), |IS) o |SI), |0T) ^ |T0), |IT) ^ |TI) 

(2) 

where the notation o denotes a mapping under ^/i^' j^x' 

for instance, /i|'j^^|OS) — 2|S0), as well as h[i^^\SO) = 
2 1 OS). Terms in Eq.d^]) are called "hopping" because 
they exchange the two rung states involved. The scat- 
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tering terms, on the other hand, change the rung states 
and are, 

Mlm(|TT)-|SS))=4(|0I) + |I0)), (3) 

/.ll,+i(|0I) + |I0)) = 4(|TT)-|SS)). 

All other terms not listed in Eqs. ([2|) and ([3]) map to zero. 
Explicitly, 

/^M+ilOO) = 0, = 0, /ill,+i|ST) = 0, 

hl^,\TS)=0, /zl|,+i(|OI>-|IO))=0, (4) 

as well as /i|^+i(|TT) + |SS)) = 0. We note that ah di- 
rect products of |T), |S), |0) and |I) are eigenstates of the 
whole rung Hamiltonian X^i^i particular, all two- 

rung states annihilated by the leg Hamiltonian, i.e., those 
given in Eq. dU, are eigenstates of hj- + /i^^ (while for 
finite J the |TT) + |SS) is not). Therefore, if one has a 
product state in which all nearest-neighbor terms are an- 
nihilated by the leg Hamiltonian Eq. or correspond 
to hopping Eq. ([2]), H will preserve the number of each 
of the 4 letters. Such states therefore form an invariant 
subspace that can be labeled by a number of each of the 
4 letters that the states contain. 

Let us list some of the simpler invariant subspaces. 
The simplest is a product state with an interchanging 
|S) and |T) states, i.e., either |STST...) or |TSTS...). 
These two states, called an S — T background, are eigen- 
states of the XX ladder, Eq. ([T]), with the eigenenergy 
2{Nt - Ns)J - AJL, where 7Vt,s is the number of the 
respective letters. From those two states one can imme- 
diately get a series of invariant subspaces by inserting 
into an S — T "background" a number of O or I states. 
For instance, L states of the form jSTIST . . .), where one 
I is at any of L possible positions, form an invariant sub- 
space. Using the hopping property of the leg Hamilto- 
nian ([2]) we can see that when H acts on such a state 
it will either move the I to one of the two neighboring 
sites, or leave it at its original place. A completely anal- 
ogous thing happens with two (or more) inserted Is in 
the S — T background. The dimension of such an invari- 
ant subspace is (^J , where A^i is the number of Is, while 
the total magnetization is Z = —2Ni. We can see that 
within one such subspace, characterized by a fixed num- 
ber of Is, the dynamics is of the same hopping kind as 
in a single XX chain. It is therefore immediately clear 
that the transport of magnetization is ballistic within 
such an invariant subspace. The speed of a propagation 
front can also be immediately read out and is u = 4 in 
the units chosen. Using |0) instead of |I) in the above 
construction results in invariant ballistic subspaces with 
magnetization Z = +2No, where A^o is the number of 
inserted Os. The above construction gave us ballistic 
subspaces with any nonzero magnetization. There exists 
though also a ballistic subspace with zero total magne- 
tization. To see this we note that if we insert a single 



two-rung state (|0I) — |I0)) into an S — T background, 
the H acting on such a state will move I and O to the left 
and to the right, one site after each application, while it 
will preserve the number of each of the four rung letters. 
Such states, which are a product of an S — T background 
and {\Ol)j^k — |IO)j,fe) at the j-th and k-th rung therefore 
form an invariant subspace of dimension (2) [20l |. again 
with ballistic dynamics, but having total magnetization 
0. There are other, more complicated invariant subspaces 
that we shall not consider in the present work [2lj . Note 
that the total dimension of invariant subspaces discussed 
in the present work grows at least as ^ 2^ (their to- 
tal relative size though goes to zero as ^ 2^/2^^ when 
L 00). 

Considering a plethora of invariant subspaces one 
might be tempted to think that the XX ladder is an 
integrable system. While the limits J = and the 
Hubbard limit are indeed integrable, for general J and 
A the model is not integrable. It in fact displays a 
standard feature of chaotic quantum systems, namely a 
random-matrix-like repulsion between nearest-neighbor 
eigenlevels [23|. Calculating eigenenergies of the XX 




1 1.5 2 2.5 3 

FIG. 1. Level spacing statistics in the XX ladder (bars) agrees 
well with the Wigner's surmise (full curve), being a standard 
signature of quantum chaotic systems. Parameters are J = 1, 
A = 0, L = 10 (similar statistics is obtained also for A — 1.5) 
and the subspace with Z = 12 and eigenvalues in the energy 
interval E £ [1,9] are used. 

ladder and removing eigenenergies corresponding to the 
above-mentioned symmetries, the nearest-neighbor level 
spacing is well described (Fig. [IJ by the Wigner's sur- 
mise, p{s) = ^sexp (— s^7r/4), that represents a good 
approximation to level spacing statistics in systems with 
orthogonal symmetry (the so-called Gaussian orthogonal 
ensemble) j22l] . 

We have therefore seen that in the XX ladder, even- 
though it is in general chaotic, there is a ballistic sub- 
space having any total magnetization. Ballistic trans- 
port, usually associated with local conserved quantities, 
has been so far observed only in integrable systems. Here, 
however, the invariant subspace is not associated with 
any local conserved quantity [23|. All our statements on 
the existence of invariant ballistic subspaces hold for any 
value of the rung anisotropy A as well as for any J (in 
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fact, one can also allow for an inhomogeneous rung as 
well as leg coupling strength). As a simple consequence, 
there exist ballistic subsectors also in the one-dimensional 
Hubbard model at any filling, including at half-filling for 
which diffusive transport has been observed ^] (there 
is no contradiction because, as we shall see, the actual 
observed transport depends on details of the coupling 
with a bath) . A heuristic picture of how ballistic trans- 
port comes about is quite simple: provided we have the 
right background state that does not interact with the 
Is - we can consider Is as being elementary excitations 
- the Is move around freely. The same holds if only Os 
are present. However, as soon as we leave the invariant 
subspace, by allowing both Is and Os, an I can collide 
with an O, resulting in a "scattering" event that modi- 
fies the background, effectively introducing other types of 
excitation (a neighboring T — T or an S — S). An inter- 
esting question is, what is the nature of transport in the 
subspace orthogonal to all ballistic ones? While we are 
not able to analytically answer this question we shall use 
large-scale numerical simulations to demonstrate that it 
is diffusive. 

There are two ways to assess transport properties, (i) 
Coupling the system with reservoirs and studying a gen- 
uine nonequilibrium steady state, by, for instance, mea- 
suring the stationary current through the system. One 
possibility is to implement the Lindblad master equa- 
tion with cold atoms. Essential ingredients have already 
been experimentally demonstrated [25,] . (ii) Another pos- 
sibility is without any external coupling to reservoirs, by 
simply preparing a nonequilibrium initial state, and then 
observing how the initial disturbance spreads. Such ex- 
periments have been performed llSl . Il4 | as well as ap- 
propriate numerical simulations '29] . For our ladder sys- 
tem one could prepare initial states from the ballistic 
subspace by quenching a strong dimerized rung coupling 
strength, so that the ground state would be |STS...), 
as well as strong local magnetic field introducing I or O 
"excitations" . 

Numerical demonstration.- Here we focus on numer- 
ically demonstrating (we set J = I) the above findings 
via a master equation of the Lindblad form 27| , in which 
the action of reservoirs is described by a dissipator C^^^, 



dp/dt = i[p,H]+C'''^p) = C{p). 



(5) 



The dissipator C^^^ is expressed in terms of Lindblad op- 
erators L^, ^'^'^(p) = Efc ([^fcP, 4] + I^fe: p4]) ■ Af- 
ter long time 2^] the solution of the Lindblad master 
equation ([5|) converges to a nonequilibrium stationary 
state (NESS) poo — hmt^oop(i). In all cases studied 
here the NESS is unique and therefore independent of 
the initial state. We are, in particular, interested in the 
NESS expectation value of the magnetization current J 
in the upper leg, J = tr(poojl'^^), jj'^^ 
^k^k+l 



for the lower leg. Crucial question then is how the current 
J scales with L, provided we keep the driving constant. 
If one has J ^ the transport is diffusive, if J ~ L", 
it is ballistic. We shall demonstrate that in the XX lad- 
der, depending on the choice of Lindblad operators, one 
can have either ballistic or diffusive behavior. 

Ballistic case.- Here we want to demonstrate ballis- 
tic transport in one of the invariant subspaces. We 
choose the subspace of dimension 2(L -I- 1) spanned by 
two background states and 2L states containing one I 
in an S — T background. The idea is to use such Lind- 
blad operators that inject I at the left end and remove 
it at the right end, while at the same time preserving 
the invariant subspace. One possible choice of Lind- 
blad operators is Li = |ITS . . . ST)(STS . . . ST|, L2 = 
|IST...TS)(TST...TS|, L3 = |STS . . . ST)(STS . . . SI| 
and Li = |TST...TS)(TS...TI| (written for even L). 
Li,2 inject one I while ^3^4 absorb one I. They therefore 
preserve a union of zero and one I-excitation subspaces. 
It should be clear from our explicit exposition in the first 
part of the Letter that within this subspace the dynamics 
should be ballistic. We have diagonalized the superop- 
erator corresponding to the Lindblad equation ([5]) with 
the chosen four Lindblad operators. Finding NESS poo 
for L up to 30, we can see in Fig. [2] that the current J 
is indeed independent of the system size. The magneti- 
zation profile is, as expected for a ballistic system, flat 
in the bulk (magnetization is the same at 1st L — 1 sites 
and different at the last; data not shown). In fact, a very 
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'^('^k'^k+i ~ '^^'^k+i)- A similar expression can be defined 



FIG. 2. Ballistic transport in the zero and single I-excitation 
subspace of the XX ladder, A = 0. Scaled current J/{zi— zl) 
is independent of the system length L. 

simple formula for the current holds, J — ^ (zi — z^), 
where w = 4 is the speed given by the bandwidth due to 
the hopping ([2|) and Zk = ti{pooCrl). Because of the sym- 
metric driving currents and magnetization are the same 
in both legs. 

Local Lindblad operators.- Here we will take eight 
Lindblad operators that will induce a nonequilibrium 
situation. They are Li,2 = -^/I ± piacf at the left 
end, L'iA — y/1 ± pLLaC^ at the right end, and, simi- 
larly, £5,6 = VI ± MirTi , L7,8 = yr±7IE7T^. Driv- 
ing parameters p, try to induce a nonzero magnetiza- 
tion at the respective ladder site. Lk do not preserve 
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any of the invariant ballistic subspaces. We shall con- 
sider two cases: (a) symmetric driving around z — 0, 
-Mia = -MiT = ^J■La = ^J■LT = and (b) an asym- 
metric driving inducing a nonzero total magnetization, 
/iio- = 0.5, ^Lcr = 0.9, /iiT = OA,fj,LT — 0.8. In the 
case (a) all contributions to the current from invariant 
ballistic subspaces mutually cancel and there is no net 
ballistic current (if there is a contribution from a sub- 
space 5, there is also one from a spin-flipped counterpart 
F(S), in which all Os are replaced by Is, and which car- 
ries an opposite current). In the case (b) there is no such 
cancellation and there is a nonzero net ballistic contribu- 
tion. For finite L, one will therefore have a combination 
of ballistic and diffusive contributions, with the ballis- 
tic one dominating long-time behavior. In (a) the NESS 
is, for small driving fi, close to the identity matrix (for 
/i = the NESS is poo oc 1). For the asymmetric Lind- 
blad operators, and if pia- ^ fJ^La — Pir = P-Lt = P-, the 
NESS is Poo cx Uf^iil + Mct|)(1 + flrj). There- 
fore, comparing NESS states with the grandcanonical 
one, p cx exp {—[i[H — (t)Z)), we can say that the NESS is 
for small driving in (a) close to the infinite temperature 
state, while in (b) it is close to an infinite temperature 
but finite (3(f> state. Different transport behaviors seen, 
e.g., in Fig. 0] are therefore exhibited close to thermal 
states. There is in principle no difficulty in using the 
more complicated Lindblad operators leading to states 
at finite temperatures (see, e.g., 2nd Ref. under Q for 
an example) for which similar results are expected. For 
small driving NESS is therefore close to separable in the 
operator space, and one of the best methods for solv- 
ing the Lindblad equation is a time-dependent density 
matrix renormalization group method in the operator 
space. Details of our implementation can be found in 
Ref. 30] . We managed to calculate NESS for systems of 
up to L = 100 sites. In Fig. |3]we show a magnetization 
profile along one of the legs (it is the same in both) for 
the symmetric driving. We note that the linear profile 
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FIG. 3. Magnetization profile along the ladder in the NESS 
for the symmetric driving (a) with 8 Lindblad operators. L — 
64, 11 = 0.1, A = 0. 



fusive J w ^(zi — zl) relation. On the other hand, for 
the asymmetric driving the scaling of the current with 
the system size is markedly different - the ballistic con- 
tribution is nonzero and dominant, circles in Fig. |4l An 
important question is, how stable are the phenomena pre- 
sented; for instance, can one still have ballistic transport 
after an addition of small perturbation to H given by 
Eq. ([T])? While a detailed answer goes beyond the present 
work, let us just mention that small interaction along 
legs of the form Jz{afa^_^_i +t^t^^i) has very little infiu- 
ence on the current. Concretely, choosing = 0.1 and 
the asymmetric driving current changes by about ~ 5 % 
and would be almost indistinguishable from the data for 
Jz = (circles in Fig. |4]). 
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FIG. 4. Scaling of the current with the system size L for 
the driving with 8 Lindblad operators. Squares are for the 
symmetric driving with p = 0.2 and A = 0, circles for the 
asymmetric one with A = 1.5 (in this case we plot J/{pl ~ 
pi) = J/0.4). The full line suggests diffusive ~ l/L scaling. 

Conclusion. - We have identified a number of ballistic 
invariant subspaces in a class of XX spin ladders. Such 
systems are in general quantum chaotic, with integrable 
limits, one of them being the one-dimensional Hubbard 
model, appearing at special parameter values. We there- 
fore explicitly show that one can have ballistic transport 
in a chaotic system. Using extensive numerical simula- 
tions we in addition show that, outside of these invari- 
ant subspaces, the transport is diffusive. In a simple 
strongly correlated system one can therefore observe ei- 
ther ballistic or diffusive transport. Such behavior can be 
observed either in the time-evolution of inhomogeneous 
initial states, or in an open-system with external driving. 
Which of the two regimes is exhibited depends on the 
initial state or, in a master equation setting, on the sym- 
metry of a particular driving. Support by the research 
Program P 1-0044 is acknowledged. 



is typical for diffusive systems. That the transport is in- 
deed diffusive is confirmed by scaling of the current J 
with the system size (Fig. 21 squares), showing a nice dif- 
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